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SINGULAR LIMITS IN HIGHER ORDER LIOVILLE-TYPE EQUATIONS 


FABRIZIO MORLANDO 


Abstract. In this paper we consider the following higher order boundary value problem 

f (-A)'"ti = p^'"y(x)e“ in a 
1^ BjU — 0, \j\ < m — 1 on 

where Q is a smooth bounded domain in with m G N, V{x) 0 is a smooth function positive somewhere in H 

and p is a positive small parameter. Here, the operator Bj stands for either Navier or Dirichlet boundary conditions. 
We find sufficient conditions under which, as p approaches 0, there exists an explicit class of solutions which admit 
a concentration behavior with a prescribed bubble profile around some given fc-points in Q, for any given integer k. 
These are the so-called singular limits. The candidate fc-points of concentration must be critical points of a suitable 
finite dimensional functional explicitly defined in terms of the potential V and the higher order Green’s function with 
respect to the imposed boundary conditions. 


1 . Introduction 


Let n C be a smooth bounded domain with m G N. This paper deals with the existence, the qualitative 

properties and the asymptotic behavior of nontrivial solutions u : SI —>■ R to the following higher order boundary value 
problem 


r (-A)"*!! = in SI 

BjU = 0, IjI < m — 1 on i9SI, 


where we prescribe the boundary conditions to be either Navier BjU := (—A)-’m or Dirichlet BjU := dlu with u = v{x) 
the unit outer normal, V (x) yf 0 a given smooth potential and p G R^ a small parameter which tends to zero from above. 
The leading part (—A)"* is, usually, called polyharmonic operator on R^™', it is simply obtained as the composition 
m-times of —A. In the sequel, we denote by Hm(SI) the linear space of polyharmonic functions of order m in SI. 
We recall that CH) corresponds to a standard case of uniform singular convergence, in the sense that the associated 
nonlinear coefficient p^"^k(x) goes to zero uniformly in SI as p 0. In the last decades a lot of work has been done in 
the study of the asymptotic behavior of solutions in the limit, the so-called singular limits. 

Higher order elliptic equations involving exponential nonlinearitles appear naturally in conformal geometry. On a 
compact surface (M, g) with a Riemannian metric g, a natural curvature invariant associated with the Laplace-Beltrami 
operator is the Gaussian curvature K — Kg. Under the conformal change of metric 'g = we have 

-Aj + K = KuC^^ on M 


where Ku denotes the Gaussian curvature of {M,g). For compact manifolds of general dimension N, when N is even, 
the existence of a N-th order operator conformally covariant of bi-degree (0, N) was verified in [52] • This operator 
is usually referred to as the GJMS operator and its leading term is the polyharmonic operator. However, this operator 
is only explicitly known on the standard Euclidean space R'^ where it reduces to its leading term and on the standard 
sphere through the Beckner-Branson formula isi nni- Attached to this operator, there is a natural concept of 
curvature, namely the -curvature, firstly introduced by Branson in [9]. If A = 2m is even, we have after a conformal 
change of metric g = e^^g that 


pN ^ -JVupiV 


( 1 . 2 ) 


and 


nN , r'.N f-.N Nu 
P„u + Qg =Qg e 


g I -^9 '^9 

We notice that (II. .Ill is a generalized version of Gauss’s identity in dimension 2. 


(1.3) 


One can think the GJMS-operator as a higher order analog of the Laplace-Beltrami operator defined on a compact 
2-manifold, and the Q-curvature can be thought of as a higher order analog of the Gaussian curvature Kg, in fact, 
in dimension 2 we simply have = — Ag and Qg = Kg. We refer to the research monograph m for related topics 
and recent developments. Therefore, in the spirit of the classical Poincare uniformization theorem for closed surfaces, 
a similar question arises: for every closed N-dimensional manifold {M, g) with N even, does there exist a metric g, 
conformally equivalent to g, for which the corresponding Q-curvature Qf is constant? 

In view of the transformation law g = e^^g, the problem reduces to solve the following A-th order nonlinear elliptic 
equation: 

g u + Qg = Qe in M, 


1 


(1.4) 
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where Q is a real constant. The above equation has been extensively studied during the last decades. In (sa. Theorem 
1.1), generalizing a result in a seminal paper by Djadli and Malchiodi in the fourth dimensional case |24| . the author 
solves completely this problem under generic assumptions: 

kerP^ = {constant}, / Q^dvg 7^ k{N — l)!a;iv for k G N*. 

J M 

Here we are interested in the special case where the manifold is the Euclidean space endowed with the Euclidean 
metric g^ 2 m. It N = 2m, Pg^„, are pointwise operators and by stereographic projection on R^"*, they simply become 
(—A)"*, whereas, if N is odd, P^;^ are nonlocal pseudo differential operators and by stereographic projection on R^, 

they simply become (—A)"2'. In both cases, by the flatness of R^, we have = 0. In the case N even, (11.411 reduces 

to 

i-A^u = in R^'", (1.5) 

which gives rise to a conformal metric g = e^^g^ 2 m whose Q-curvature is given by the constant Q. Clearly we can set 
(lESl) in a regular bounded domain Q, C R^"* but, since we are no more dealing with a manifold without a boundary, 
we need to impose some boundary conditions on dQ. in order to have a well-posed problem. The main analytic difficult 
in working with nonlinear higher order elliptic equations lies, basically, in the failure of Maximum Principle and all the 
comparison type results based on it. Moreover, it is well known that for higher order problems the lack of the Maximum 
Principle strongly depends on the kind of boundary conditions imposed on the solutions. In the case of Navier boundary 
conditions the positivity property is not sensitive to the geometric or topological characteristics of the domain. On the 
other hand, if we consider Dirichlet boundary conditions, this property is guaranteed only in some special domains, for 
example in balls taking into account the classical Baggio’s formula [7]. Recent results show that it does not hold for 
certain ellipses |30l I54| and for squares |2n| but it holds for domains close (in a suitable sense) to the planar disk [33], 
or some non-convex domains m- We refer to the research monograph m for further details and historical information. 

The general behavior of arbitrary families of blowing-up solutions to problem (El) when infn V > 0 has become 
understood after ([44|, Theorem 1) in the case of Dirichlet boundary condition and after ([38], Theorem 1.2) in the case 
of Navier boundary conditions. It is known that if Up is an unbounded family of solutions for which V(x)e“'’dx 

remains uniformly bounded, then necessarily 


1* ^ 
lim p 

P-^O 


^ f Vix) 

Jn 


e^^dx = kA2r} 


( 1 . 6 ) 


for some integer A: > 1, where A2m := 2^"‘m!(m — l)!tJ2m is a normalization constant and uJ 2 m denotes the area of 
the (2m — l)-dimensional unit sphere in R^"". Notice that the constant A2m has a geometric meaning: it is the total 
Q-curvature of the round 2m-dimensional sphere, see also for further information. Moreover, there are fc-tuples 

of distinct points of Q, the so-called concentration points, {^ 1 , separated at uniformly positive distance from 

each other and from the boundary dtl as p }, 0 for which Up remains uniformly bounded on D \ U2=i^'5(G') ^^d 
sups^(5.) Up —>■ -1-00 for any 5 > 0. 

An obvious question is the reciprocal, namely construction of solutions to Problem El with the property (II.6II . 
We shall see, firstly, that the location of the concentration points is related to the set of critical points of the finite 
dimensional functional tpk defined explicitly in terms of the potential V and the Green’s function G of the polyharmonic 
operator with respect to the Dirichlet or Navier boundary condition respectively by 


=: -E [21ogV(0) +7I(G',C.)] (1-7) 

defined for points ^ = (Cij £,k) £ fl*’ \ A, where fl*’ denote the cartesian product of k copies of and A := {^7 £ fl*’ : 
Ci = G', * 7^ J} denotes the diagonal in fl*’. For any ^ G fl, let G—G{x,^) denote the Green’s function of (—A)™' under 
the appropriate boundary conditions as the unique solution to 

(-ArG(;0 = A2mSdx) infl 
BjG{-, G = 0, \ j\ < m — 1 on dO,, 


(1.8) 


where S^{x) := S{x — ^) is the Dirac measure centered at the pole G Clearly, the concept of Green’s function generalizes 
that of a fundamental solution. If the differential operator involved has constant coefficients, as in our case, it is in 
fact often advantageous to think of the Green’s function as a perturbation of the fundamental solution. Namely, we 
decompose G in a singular and regular part as 

G{x,0 = Hix,0 + K{x,0 (1.9) 


where K{x,ff) = 4m log is the fundamental solution of (—A)™' in R^’ 
auxiliary polyharmonic function H{x,^) so that 

(_A)'"H(-,G =0 infi 

ondn. 


and the regular part of G is a smooth 


( 1 . 10 ) 
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Finally, let H{^) denote the diagonal of the regular part H, usually referred as the Robin’s function of Q at 

^ and, notice that —oo as ^ dfl. 

We observe that arbitrary critical points of ipk are not all candidates to be concentration points. Since the construction 
relies on perturbation technique those critical points must also satisfy a sort of nondegeneracy condition. In [3, the 
authors solve the question in the case m = 2, V = 1 and with Navier boundary conditions. They established that for 
any nondegenerate critical point of ipk, a family of solutions concentrating at this point as p 0 does exist. But, as 
remarked in [^, their construction, based on a very precise approximation of the actual solution and an application 
of Banach fixed point theorem, uses nondegeneracy in essential way. This assumption, however, is hard to check in 
practice and in general not true, an annulus being an obvious example. Next, in m the authors present a construction 
of a blowing-up families of solutions under a weaker nondegeneracy assumption of [5], namely that ipfe has a topologically 
nontrivial critical value and including also a general smooth potential V. 

In the same spirit of m , we provide sufficient conditions for the existence of multipeak solutions to (El), generalizing 
all the previous results to the polyharmonic case by also including the Dirichlet boundary condition. We emphasize that 
all the results presented so far in literature based on the hnite dimensional reduction technique are all related to the 
Navier boundary conditions, see for instance 1261 l34l 142114^ . The critical points we will deal with are those that can 
be captured in a general way with a local min-max characterization. We consider the role of non-trivial critical values 
in an appropriate subset of a functional ip in existence of blowing-up solutions to El in order to relax the stronger and 
hard-to-check nondegeneracy condition. This local notion of non-trivial critical value was hrstly introduced in |22| in 
the analysis of concentration phenomena in nonlinear Schrodinger equation and it is widely used since. More precisely, 
we consider the following setting. A hrst observation we make is that in any compact subset of , we may define, 
without ambiguity, 

iptei = -00 if for some i 7 ^ j. 

We shall assume that there exists an open subset T> of with smooth boundary compactly contained in and such 
that infu 1/ > 0. Let ip : (D —>■ R be a smooth functional. We will say that tp links in T> at critical level C relative to B 
and Bo if B and Bq are closed subsets of B with B connected and Bo d B such that the following condition holds. 

Let us set T to be the class of all maps >I> G C[B, T>) such that there exists an homotopy ^ G C([0,1] x B, B) satisfying: 

^(0, •) = ids, <11(1, •) = $, ^(t, ■)\Bo = idBo for all t G [0,1]. 

We assume additionally 

sup p{y) <C= inf sup (^($( 1 /)), (1.11) 

yeBo y^B 

and for every y G dB such that p{y) = C, there exists a direction vector Ty tangent to dB at y such that a transversality 
condition holds: 

dTyP{y) = '’^p{y)-Ty^Q, ( 1 . 12 ) 

where d-ry denotes tangential derivative. We observe that this condition can also be stated in the level set framework 
as {y G dB> : p{y) = C} = 0. Furthermore, we call the min-max value C a non-trivial critical level or a stable critical 
level of the functional p in B. We observe that in the standard language of calculus of variations, the sets B, Bo link 
in B. Note that under these assumptions the min-max value C is a critical value (inside fl) for p in B, which is in some 
sense topologically non-trivial. Thus these conditions ensure the existence of a critical point of p at level C, i.e. it exists 
y d B critical point of p at level C, that is with p{y) — C. In fact, the condition ( 11 . 111 ) is necessary in order to ’’seal” 
B at level C so that, exploiting the local linking structure described, by standard deformation argument involving the 
negative gradient flow of p we are able to provide the presence of such critical point at level C in B, possibly admitting 
fully degeneracy. Not only this, any function ’’C^-close” to p (i.e. a small perturbation of p) inherits such a critical 
point. Condition (11.111) is a general way of describing a change of topology in the level sets {p < c} in H taking place 
at c = C. The min-max value C is a critical value if B is invariant under the negative gradient flow of the functional, 
while if this is not the case, one can use condition ( 11 . 111 ) to prevent intersection of the level set C with the boundary 
in order to modify the gradient vector field of p near the boundary of B at the level C and thus obtain a new vector 
field with the same stationary points and such that B is invariant and the functional ip is a Lyapunov function for the 
associated negative flow near the level C. As an example, taking B — Bo = dB, it is not hard to check that the above 
conditions hold if 

inf p{x) < inf p{x), or sup (p(a;) > sup p{x), 

xeT> xed-D j-gi, x&a-D 

namely the case of (possibly degenerate) general local minimum, maximum points or saddle-points of p. The level C 
may be taken in these cases respectively as that of the minimum and the maximum of p in B. These holds also if 
p is ”C^-close” to a function with a non-degenerate critical point in B. We will show that, for every A: < 1, the set 
B ~ {x d Q ■. dist{x, dfl) > 5} has the property (II.fill at a given C, for S small enough. This allows us to prove our 
main result. 
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Theorem 1.1. 

Let k > 1 and assume that there exists an open subset T> of with smooth boundary, eompactly contained in , 
with infc y > 0 where (fk has a non-trivial critical level C. Then, for p small enough, there exists a solution Up to (HU 
with 

lim [ V{x)e'^'’ = fcA 2 m. 

Jn 

Moreover, there is an k-tuple (^i, (,k) G T), such that as p —>■ 0 

Vipfc(^i, ^k) —>■ 0, Cfe) —>■ C, 

for which Up remains uniformly bounded on Ll \ U^=i and, for any 5 > Q, supaj( 5 ^.) Up +oo. 

Notice that for the two-dimensional version of problem (EH) many authors, using different perturbation techniques, 
have constructed solutions which admit a concentration behavior with a prescribed bubble profile as indicated by the 
blow-up analysis in the seminal papers fm 15511551 HI] . In a, providing that is not simply connected, the authors 
showed that for any non-degenerate critical point of the reduced function a sequence Up of solutions is constructed, 
which converges to a function u* = in \ S). After, in |23[ I27| the authors, independently, 

generalize this result by relaxing the assumption of non degeneracy into a weaker stability assumption, the topologically 
nontrivial critical value used also in our paper. 

The paper is organized as follows. Section [2] is devoted to describing a first approximation for the solution and to 
estimating the error. The predicted solutions are found as a small perturbation of this initial approximation. Fur¬ 
thermore, problem (11.11) is written as a fixed point problem, involving a linear operator. In Section [3] we study the 
bounded invertibility of the linearized operator in a suitable L°°-weighted space. In Section [4] we solve a projected 
non-linear problem. In this Section, in order to include Dirichlet boundary conditions we are able to avoid the use 
of Maximum Principle in the proof of Lemma This result represents the major technical difference with respect 
to m- In Sections [5] and [6] we show that solving the entire non-linear problem reduces to finding critical points of a 
certain functional and this correspond to adjust variationally the location and the high of the bubbles. Finally, Section 
[3 is devoted to the proof of the main result. 


2. A FIRST APPROXIMATION OF SOLUTION 


In this section we construct a reasonably good approximation U for a solution of EH, we give an expansion of its 
associated energy and then we estimate the error of such approximation in appropriate norms. Our construction relies 
in an essential way in two main steps: 

• on the choice of a suitable family of approximating solution, 

• on the invertibility, in some sense, of the linearized operator evaluated at such approximating solution. 

A useful tool involved in this construction is, in fact, concerned with the classification of entire solutions to the higher 
order Liouville equation which represents, essentially, the associate limit problem to EH- Given a constant Q G R, we 
consider the following limit profile problem 


(-A)’"-a= (36=™-“ 
u{x) = o(|a;p) 


in R^^-.m > 1 
as |a;| —>■ oo. 


( 2 . 1 ) 


In [45] it has been proved that all the solutions of EH are radially symmetric about some point ^ G R^™' and are all 
given by the one-parameter family of functions 


UsA^) 


log 


((52 -G I® - 


( 2 . 2 ) 


for any free parameter (5 > 0, ^ G R^"* and where a 2 m ■= 2^"‘"'"^m(m — 1)! if and only if the constant Q is strictly 
positive as pointed out in [3B] , if Q < 0 and m > 1 there are no standard solutions and there are no solutions presenting 
a nice behavior at infinity. We call the functions of this form standard (or geometric) solutions. Notice that given 
a solution u to {—A)"'u = Qe^""^ in R^™',m > 1 and A > 0, the new function defined as v := u — ^logA solves 
(—A)"*u = XQe^""^ in R^"*, hence what really matters is just the sign of Q. There is no loss of generality if one 

assumes that Q G {0, ±(2m — 1)!}. On the other hand, the particular positive choice Q = (2m — 1)! has the advantage 
of being exactly the Q-curvature of the round sphere 5^"“. This fact implies that the standard solutions has the 
geometric property that = (7r“^)*gg2m, where tt : > R^"* is the stereographic projection, then by Mobius 

transformations actually give us the large family of solutions EH- We recall that in dimension 2, by employing the 
method of moving plane, in [18] the authors were able to classify all the radial solutions of 

r (-A)m = 6^“ in R^ 

\ e^“ G L^(R^) 


(2.3) 
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as standard solutions without any other condition. Conversely, if m > 1 a careful study of radial solutions shows that 
there are solutions which do not come from the smooth function on 5'^"* through the stereographic projection. Thus, to 
overcome this difficulty, first in |58| the author added a constraint in the behavior at infinity of the solution in order to 
have the same classification obtained in m = 1 case. Notice that their proof appears to be overly simplified, for instance 
their Lemma 2.2 is not conclusive. In f|45|. Theorem 2) the author give a proof of this result with all the details and 
develop some criteria to characterize of the non-standard solutions. Of particular importance is the following criterium: 
If a solution u is non-standard, then there exist 1 < j < m — 1 and a constant a ^ 0 such that 

lim A^u(x) = a. 

\x\—^oo 

Notice that a similar property can be shown for every solution in the negative case. Moreover, the standard solutions 
satisfy some nice properties: 

(1) A2m {mass quantization), 

(2) A2m5c in measure sense as 5 —>■ 0 {concentration property), 

(3) given a small (5o > 0, sup^^ Us,^{x) —>■ -boo as 5 —>■ 0 and Us,^ remains uniformly bounded on \ Bsq{^) 
{Us,^{x) blows up at ^ and remains bounded away from it), 

(4) Us,i{x) —>■ —00 as |a:| —>■ -boo. 


Due to these all properties we shall use Us,s,, suitable scaled and projected, as fundamental building blocks to construct 
an approximate solution of dJ around Given k a positive integer, let us consider k distinct and well-separated 
points ^ 1 , ...,^fe G where the spikes are meant to take place with V{^i) > 0 for all i = 1, ..., k because, as we said 
before, for Q < 0 and m > 1 there are no standard solutions at all. As we will see, a convenient set to select these 
points is 

Af G : dist(^i, 9D) > 2<5o Vi = 1,..., k, min |^i — Gl ^ 2(5o} (2.4) 

where 5o > 0 is a small but fixed number. From now on we fix ^ = (^i, ■■■,^k) G M. Such a definition is motivated by 
a blow-up analysis performed in. Given e > 0, we define for x G around each 

(l + eb 


SO that Ue is a rotationally symmetric solution to 

(-A)"*u = p^’"l/(Ci)e“ in 

with 

2m _ a2m(2m-l)!e2- 


(2.5) 


( 2 . 6 ) 


u (l + e2)2- ’ 

that is, p ~ e as e 0. Now, let us notice that equation ( 1 ^ is scale invariant under some dilation in the following 
sense: if Ug is a solution of ( 1 ^ and Pi > 0 , i = 1 ,..., k, are k positive parameters to be selected properly later on, 
then We(pi-) -b 2mlogPi ia also a solution of (12.61) . Thus, with this observation in mind, we define for all pi > 0 

Pi (1 -b e^) 


SO that the solution Ug will look, near each ^i, like Ui for certain e- independent parameters pi. 

Remark 2.1. 


(2.8) 


(1) Note also that Ui uniquely solves (12.61) in entire R^"* with p e as e —>■ 0, in fact, up to scaling and translation 
invariance, these are the only solutions with finite energy condition and a nice behavior at infinity. 

(2) A considerable characteristic of elliptic PDFs with nonlinearity of exponential type in the right dimension (i.e. 
(—A)*” in R^"*) is that all the bubble solutions have proportional heights with coefficient of proportionality pi. 

We would like to take first approximation to a solution of (11.11) . Since, obviously, BjUi, where \j\ < m — 1 

and i = 1 ,..., k, are not zero on the boundary 9D, we perform a suitable polyharmonic correction to modify in 

order to fit in the boundary conditions. Let Hi{x) be the smooth unique solution of 

r {-A)’-Hi{x) = 0 inn 

\ BjHi{x) = BjUi{x),\j\ < m — 1 on dfl. ^ ' 


We define our first approximation U{^) as 

k 

U{^) := X] C/i, Ui := Ui + Hi (2.10) 

i=l 

which now, by definition, satisfies the boundary conditions. We want to find out the asymptotic behavior of Ui away 
from and around ^i. 

First, we recall the standard regularity statement and an a priori estimate for higher order linear problems in the U" 
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space framework developed by Agmon, Douglis and Nirenberg [T] about 50 years ago. Before this, we remark that, in 
general, the choice of the linear boundary operator Bj’s for polyharmonic problems is not completely free. In order to 
obtain a priori estimates and, in turn, existence and uniqueness results we need to impose a certain algebraic constraint, 
the so-called Shapiro-Lopatinksii complementary condition. For any j, let Bj denote the highest order part of Bj which 
is precisely of order mj. For every point x £ Q, let n{x) denote the normal unit vector. We say that the complementing 
condition holds for BjU = gj on dfl if, for any nontrivial tangential vector t{x), the polynomials in t Bj{x;T + tv) 
are linearly independent modulo the polynomial (t — ilrl)"* with i the unit imaginary number. We remark that both 
Dirichlet and Navier boundary conditions satisfy this condition. For example, the homogeneous Neumann boundary 
conditions Au = dv{Au) = 0 on dO. for the biharmonic equation A^ti = 0 in not satisfy it. 

Lemma 2.2. 

Let Q, a bounded domain in 
following problem 


for all f € for all gj £ with mj £ N the maximal order of derivatives of Bj and assume 

that Bj satisfy the Shapiro-Lopatinksii condition. Then (12.1111 admits a unique strong solution u £ Moreover, 

there exist a constant c = c(|f2|, k, m, Bj) >0 independent of f and of the gj ’s, such that the following a priori estimate 
holds 

||u||iv‘i.p(n) < ll/llw9-2m.p(n) + (2-12) 

\ .—: W V (c?nj / 

^ J=1 

Now, arguing as in the classical [^, we obtain the following crucial characterization in which we asymptotically 
expands Ui in LI. 

Lemma 2.3. 

Assume ^ £ M. Then we have, as e —^ 0 ; 

(1) uniformly in LI in the -sense 

Hi{x) = H{x,^i) - 2mlog/q(l + e^) -|- logy(^i) + 0{ple^), (2.13) 

(2) uniformly in the region \x — ^i| > (5o 

Ui{x) = 2mlog/ii(l + e^) - log V'(fi) - 4m log I* - ^i| + (2.14) 

(3) In particular, in this region, as we expect by the asymptotic analysis, 

Ui{x) = G{x,^i) + OifiU^) (2.15) 

where a is an arbitrary Holder exponent, G{x,^) is the Green’s function and H{x,^) its regular part. 

Proof 

Let us prove (12.131) . Define z{x) := Hi{x) — H{x,f,i) -\- 2mlog/ri(l -I- e^) — log V(^i), then z £ HmiLl) and it satisfies 

(-A)"*«(a;) = 0 in D 

z{x) = —Ui{x) — 4m log I ■ —^i| + 2m log Pi (1 -I- e^) — log F(Ci) on dLl 
Bjz{x) = BjUi[x), |j| < m — 1 on dLl. 

Since Bjz{x) = 0{pL^e^) with \j\ < m — 1, using the standard L^-estimates in Lemma 12.21 and the Sobolev embedding 
theorem we get z{x) = 0{g,\e^) uniformly in in the Gifff~ ’“-sense. 

The point (12.1411 follows, easily, from the definition of Ui in (12.811 . Finally, combining (12.1311 and (12.141) we get (12.1511 . 

□ 

While Ui{x) is a good approximation to a solution of HU near ^i, U{^) is a good approximation far from these points 
but, unfortunately, it is not good enough for our constructions close to these points. This is the reason why we need to 
further adjust our ansatz. To do this, we need to improve that the remainder term U — Ui(x) = (Hi + Uk) ~ 0, 

that is it vanishes at main order near as e —>■ 0. Using Lemma [2.3l we can achieve this through the following precise 
selection of the parameters pi 

2mlog Pi = iy(6,6)+logU(^i) + E^(^Tei)- (2.16) 

We thus fix Pi a priori as a function of ^ £ Ai, fXi = tJ-i{(,) Vi = 1,..., k. Since ^ £ M, for some constant c > 0, 

i < Pi < c for all i = 1,..., k. (2.17) 

c 

and, thus, instead of O(pie^) we can write O(e^) in all the following asymptotic expansions. 

Now, let us denote with 


Let 1 < p < oo and fix an integer q > 2m. Assume dLl £ C"^. Consider the 
\ {-^Tu{x) = f inLl 

\ Bju{x) = gj{x),\j\ < m — 1 on dLl, ^ ' 


Lie ■■= e~^Ll = {y£ R^"" ■ ey £ LI} 


(2.18) 
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the expanded domain with the scaling measure property |17el = e and 

(2.19) 

the expanded variables. A useful observation is that u satisfies (HU if and only if 'w{y) := u{ey) + 2mlogpe solves the 
equivalent problem 

r = V{ey)e'^ in ft, 

< w = 2mlogpe on dft, (2.20) 

(_ BjW — 0, |j| < m — 1 on dfl,. 


Let us define W{y) ;= U{ey) + 2m\og pe as the first approximation of (12.2011 with U our approximation solution (12.1011 . 
We want see how W behaves, namely, we want to measure how well it solves (12.2011 . It is convenient to do so in terms 
of the following L°°-weighted norm defined for all w G L°°{fl,) 


:= sup 


E 


(i + iy-air 


I ^ 

+ e 


k(j/)l- 


( 2 . 21 ) 


Thus, as anticipated, in the next lemma we measure the size of the error of approximation i?(y) := {—A)"'W — V{£y)e^, 
proving that although is not zero it is small in the sense of the norm defined above. We remark that the suitable choice 
of the parameters pi in (12.211) is done just to make the error term small. 


Lemma 2.4. 


Assume ^ G A4. Then we have for any y G ft,, as e ^ 0; 

||i?||.<C£ (2.22) 

where C denotes a generic positive constant independent of e and 
Proof 

Let us assume first \y — ^'n\ < for some index n. We have 

{-ArW{y) = {-Ar{U{ey) + 2mlogpe) = i-ArU{ey) = e^^[i-ArU]{ey) 

because the higher order Laplace operator is linear and it scales like a homogeneous polynomial of degree 2m. Now, 
recalling (12.101) and the fact that Hi G 'Hm(fle) we get 

£^"^[{-ArU]i£y) = e^"^[{-Arui]i£y) = ^ 


I 2\2m 

1 2m T T (j- \ H'n v t I 2m 2m \ ^ t /■ / «^ \ l 


ii^y) 


e^^a2m{2m-l)\ 




(1 +e2)2m g4m(^2 j2 + |y _ 


I 2m 2m \ ^ t / / /- 

+ e P 2^V{i 

i^n 




= a2m{2m - 1)! 


pl^{l+ey 


{ple^ + \y - i^r,) 


I 2m 2m \ ^ ^ \ ua (ev^ 

+ £ P 2^V{ii)e 


i^n 


,,'2m(-] I 2\2m 

= ^^rn{2m - + o{e^n 


(2.23) 


' (Pn£^ + \y- (Cn) 

because of i 7^ n means that y is away from the points and then we should use (I2T41) . Let us estimate V{ey)e^^^\ 
By (12.131) we have 

Hk{x) = - 2mlog/r„ + log’L(Cn) + O(e^) + 0{\x - (f„|) 

and by the definition of pi ’s (I2T61) we have 



= - E G(6, Cn) + 0(e") + 0{\x - e„|) 


i^n 

and, if i 7^ n, by (12.151) we have 

Ui{x) = G(^i, eO + O(e') + 0{\x - ^„|). 

Then, by adding, we have 

Hr.{x) + ^ Ui{x) = 0(e") + 0{\x - e„|) 


ij^n 

Therefore, 



V{ey)e^^^'> = V{£y){pef'^e^^^^'> = I/(ei/)(pe)^"‘exp \^„{ey) + Hr,{ey) + ^Ui{ey) 

^ i^n 


(2.24) 
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and using 


= V{ey){pef"'exp (un{£y) + O(e^) + 0{e\y - 


= V(ey)(p£) 


I 2\2m 

2m Pn [^ + £ J 


(ple^ + Ip - anv(^n) 


and, finally, 


= a2m(2m - 1)!- 


' (Pn + Ip - 

because by first order Taylor expansion we have 

V(£p) 


1 + 0{e^) + 0{ely - Cl)j 

1 + O(e^) + 0{ely — ^n|)^ 


We can conclnde that in this region 


V{^n) 

|i?(y)l < Ce 


= 1 + C>(e|i/-Cnl)- 


and thus ||_R(i/)||* = 0(e) as e ^ 0 if j/ G 

e 

Now, we assume ly — ^(1 > ^ for all i = 1 ,..., k, scaling back this is equivalent to say \x — ^i| > 5o so, using (12.1511 . 

{-^rW{y) = e^”^{-ArU{£y) = e^^{-Ar{G{ey,^i) + 0(e") = 0{£*^) 
and V{ey)e^^'^'^ = 0(e'‘’"). Hence, in this region, 

R{y) = 0(e^’") 

so that, finally, 

ll^(j/)ll* = 0(e) as e —>■ 0 in all He 

and the proof is concluded. 

□ 

Observe that problem uni) has a variational structure, in the sense that (weak) solutions for (11.11) correspond to critical 
points of the following nonlinear energy functional 

1 


Jp(u) ■- 


|(—A) 2 u\ dx — p 


I 


V(x)e^dx 


(2.25) 


where u acts on the Hilbertian Sobolev space H = {u £ H'^(Q.) such that BjU = 0 on dil, |j| < m — 1}, the natural 
settings for the boundary operator involved. In particular, in the case of Dirichlet boundary conditions we have 

H = HS^(n), 

i.e. H is defined as the completion of (7^(11) with respect to the standard Sobolev norm ||(—A)~M||^2(n) defined 
below, while in the case of Navier boundary conditions we have 

m — 1 


H = Hn '■= G H’™'(H) such that (—A)-’u = 0 on 9H, j = 0,1,..., 

H inherits the Hilbert space structure from H™'(H) when endowed with the scalar product 

(—A)^u(—A)^v dx if m = 2k 


(u,v)h ■— 


V(—A)^uV(—A)^V dx if m = 2fc + 1 


and where we use the following notation: 

(-A)^m := 


(—A)*^m if m = 2fc 
V(—A)*^m ifm = 2fc + l. 


The scalar product defined above induces the following corresponding norm on H: 

||(-A)'=M||i2(n) G R if m = 2A: 

||V(-A)'=m|U 2(„) G R""* ifm = 2A: + l. 


||u||fr := (u,v)^ = 


For general p G (l,+oo), one has the choice of taking the L^-version of this norm or the equivalent Sobolev norm 

||u||^ym,P(Q) := ||D™'u||j,p(n). 

We end this section giving an asymptotic estimate of the reduced energy Jp(U), that is the energy functional de¬ 
fined above calculated along our ansatz U. Instead of p, we use the parameter e, related to p by (irri) . to obtain the 
following result: 
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Lemma 2.5. 

With the election of fii ’s given by 

Jp{U) = loge| - Ab^k + 0(e) 

where ^ G M, ipk is the functional defined in (11.711 and bm ~ — l)!7r"*. 

Proof 

We consider the energy evaluated at U, that is 

= + |(-A)^[/i||(-A)^t/^|-p^”‘ / Vix)e 

^ Jn ^ Jn Jn 


(2.26) 




= Ii + I 2 + P- We will first evaluate the quadratic part Ii + p- 

Let i be fixed, we note that (—A)’^Ui(x) = (—A)”*(Mi + Hi) = (—A)'^Ui(x) = in 12 and BjUi = 0 on 912 

where * = 1,..., k and j = 0,..., m — 1. Then, integrating by parts, 

k 

I „ . 

h = 


I 2 = 


f e-^Ui, 

i=l 


First we expand 7i, using (12.1311 and (12.711 and recalling that Hi G T-LrniP) we have 


^ i^i 

^ (O f 

= -a2m(2m-l)!^y^ 


ttJn + 


- 2m log(pfe^ 


Now let Jo > 0 be small and fixed, independent of p. Let us define the change of variables 2 ; = + ^ney, where 

X G Bsg (^i) and y G = B_s^ (0) with the relative volume element dx = {)j,ie)^’^dy. As usual, we can split the 

above integral on 12 into two pieces, a sum of one integral on B^ (0) and one on its complement. Since the piece on 
the complement is small, precisely an 0(£^”*), it is enough expand the integral only on B_5o_ (0). 


h = -a 2 m( 2 m-l) 


fc ^ 




“7s^(0) (p?e2+p,|£y|2)2 


- 2mlog(pi£^ + |pi£2/|^) 


2 


+ H{^i + jj.iey,^i) + 0{e )'^ + 0{e ). 

Bringing out the term (pf£^)^™' and using first order Taylor expansion: 


we have 


H{^i + Pi£f/,6) = H{^i,^i) + 0{yi\ey\), 

h = ^a2m{2m - 1)\J2 ( 0 ) (1 + { ~ + \Pi£y\^) 


+ H{^i,^i) + 0{y,i\£y\) + 0{£ ) > + 0(£ ). 


By the fact that —2mlog(pi£^ + \fii£y\^) = —2mlog(pf£^)(l + |yp) we have 

- 2mlog(pf£^)(l + IpI^) 


7r = ia 2 m( 2 m-l)!^ J 


dy 


(0) (l + |yP)^ 


+ H{^,,^i) + 0{£\y\) + 0{£^) + 0{£^n- 


Now, since 


L 


dy 


(0) (l + lf/P)"'" 


0{s^) = c O(e^), 
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because, by improper integral test, the considered integral is bounded and 


L 


dy 


:0{e\y\) = f 
J B 


\y\dy 


(0) + 


_(o) (1 + _ 

Mi £ Mi £ 

because, by a comparison test, the considered integral is bounded, thus 

k 


0(e) = c 0{e) 


h = ia2m(2m - 1)! y] f 


+ +0(£)- 


dy 


(0) (l + |yP)^ 


- 2mlog(^fe^)(l + |i/|^) 


At this point, we can split the integrals involved in a piece over and one over R^"‘\i3_£Q_ (0). By direct computation, 

Mie 

2m ' 


we will show that the integrals over \ (0) are small. We have 

Mie 


/. 


dy 


and 


/ 

Jr 


lR2m\S^(0) (1 + 

iog(i + \y\‘^)dy 

K^"‘\S_£2_ (0) 


= 2muj2j 


— 2 mUJ 2 r 


r 

n I 

Jli 

M 

Jli 


-dr = 


(^»£) 


log 


So 


dr = c 


= ©(e""*) 


/ + 00 
«(L 


„2m+l ’ 


integrating by parts and recalling that 5o is small and we have 


log r 


/ 


<y^2?n"|” 1 


+ 00 «_j_oo 


/ + 0 
!o_ 


^ dr = logs). 


^2771-1-1 rp 


Mie 


Now since, clearly, e < and e < logs we have 


h = -Q2m(2m-1) 


k p 

'Si 


dy 


M2m (1 + 11 / 12 ) 2 ™- 


- 2mlog(/r?e^) - 2mlog(l + lyp) 


+ H{^i,^,)} + 0{e) = -a2m{2m-l)\ 


L 


dy 


^2m (1 + |?/P)2™ 




k p 

Am log /Tie) — 2m / 
.-—1 


log(l + |i/P)dy 


(1 + |2/P)= 


+ 0(e) 


and, hnally. 


h = -a 2 mr:'^(m - 1)! '^(H(^i,^i) - Amlogfue) - fca2m7r™'(m - 1)! + 0(e), 


where we have used the explicit values (see Appendix) 


Co = 


1. 


dy 


_ _ 7r"‘(m- 1)! 

g2m (1 + |//|2)2™ (2m-1)! 

log(l- + \yf)dy ^ 7r’”(m - 1 )! 

/K2m (1 + |2/P)2™ m(2m-l)!' 

Now, we consider the mixed term quadratic part I 2 , 


and 


Cl = 


I. 


h = -P2 


/ 

JBsAij) 


e“^C/^ = E2^EUC0[ [ 


e“^;7,-+ 




iJn\Bsg(ei)uBsg(ej) 


+ 




M,- T 7 

e ^Uj 


\p2mY.V(^i) 




/l + ^2 + I3 


Clearly, li = 0(e^™'). Since fg^ ^ = O(e^) then 


I 2 = 0(e^) 


L 




+ //(*, g) + o(£^) 


= 0(£ ) + 0 




log 




= 0(£ ) 


(2.27) 
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[ + / log(l® - 61) = 0(1). 


Now, 




p?"‘(l+e6^"‘ 


G(*.6) + 0(e"), 


by first order Taylor expansion G{x,^j) = G( 6 , 6 ) + 0{\x — 6 |) and with the same change of variables above we have 
1 


p^"* ^ k{^i)h = ia2m(2m - 1)! ^ f 

i¥=j -®_£2_ 

fJ.i£ 


Since 


then 


Since 


L 


dy 


[ 


(0) (i + |yP) 

dy 


G(6,6) + 0(e|y|)+0(£"). 


Ib^(o) (1 + lyP)"”* V™ (1 + |yP)2- 

ip2- ^ - 1)! ^ G(6.6) + o(e f 

^ ^ V Jb 


+ 0(£ ) — Co + 0(£ ) 


*7^1 


(0) (i + lyP)""* 


+ 0(s^). 


= 0(1) then O 




(0) (1 + lyP)^"* 


= 0(e). 


_^(o) (1 + lyP)^” _ 

Adding all this information we arrive to 

h = ia 2 m 7 r™'(m - 1 )!^G( 6 , 6 ) + 0 {e). 




Finally, we consider I 3 . 


73 = -p2- f V{x)e^ = -p^*" ^ f / V{x)e^ + f 

Jo. j_l BsQi^i) 


^\B5q Ui) 


V{x)e^ 


since 


V{x)e^ = 0{e^^) then 


k p 

I 3 = -p2- V / l/(a:)e“^(0+«i(-)rf3. o(e""‘). 

*=i 


Using the change of variables above we have 

k p 

6 = V-e/ 

Jb 5 (0) 


U(^i + Pi£p)exp <1 Ui (6 + k-iBy) + Hi{^i + Piey) \dy + 0(£^’"). 


Using (fUTOll we have 


/3 = -p^"*E 


(6+M» y)(^2g2 + |^.£y|2)2,ny(^.) 


Mie 

exp |77(5i + fj,iey,^i) - 2 mlogpi(l + + log U(^i) + 0 {e"‘) \dy + 0 (£"”"), 

By a first order Taylor expansion for H{^i + fj,iey,^i) we have 

U(6 +M»£1/) 


h = -p^^Y. 


(0) (pie)2-(l + |yP)2- 


exp) 1 ^( 6 , 6 ) + 0 (|pi£p|) + 0 (£^) + 0 (£^™') = 


■E 


-a2m(2m - 1)! 


(1 + £2)2 


U (6 + k^iBy) 

(0) (pi£)2™'(i + |yP)2 


exp) + 0 {e\y\) + 0 (£^) \dy + 0 {e^^) 


(2.28) 
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since (1 + ~ 1 + 0(e^™') then we have 

73 = -a 2 ^( 2 m-l)!^ [ Jjfi + + Q(g|yl))(l + 0 {s^))dy + 0{s^-^). 

Mie 

By first order Taylor expansion for V{^i+^iey), the fact that (l+0(£|i/|))(l+0(e^)) = l+0(e^)+0(£|i/|)+0(e^+e|j/|) ~ 
1 + 0 {e), the election of pi’s (12.1611 and the integral value co we have, finally, 

h =—a 2 mn’^im — ly.k + 0{e). (2.29) 

Thus, summing up (12.271) . (12.281) . (12.291) . we can conclude the following expansion of JpiU) 

p fc " 

Jp{U) = ia2m7r™'(m-1)! ^CO - 4m log/ii£ + ^ G(^j, ^i) 

i=l 


li 


— 2fcQ2m'7I'™'(m — 1)!. 


We focus on I a. 


h = -a2m7r"‘(m - 1)! 


^ 77(Ci, 6) - 4m log /ii - 4m log £ + ^ G(Cl, CO 




by the election of fii ’s (I2T61) we have 


h = -027^77"“(m - 1)! 


77(Ci, Ci) - 2 log Vm - 277(^7, Ci) - 2 ^ G(6', CO 




1 ^ 

+ -a2m7r"*(m - 1)! ^ G(Ci,CO “ 2ma2m7r"‘(m - 1)! ^ | log£| 


2 

= ia2m7r"*(m - 1)! 


k 


^-77(Ci,C0-21ogl/(e0 


ia2m7r"*(m - 1)! ^ G(Ci,CO “ 2mfca2m7r"*(m - 1)!| log£|. 


Introducing the function tpk{0 dehned as in (11.71) and the new constant bm dehned as in the statement we conclude 
the proof. 

□ 


In the following, we will stay in the expanded variable y £ Qe- 

We will look for solutions of problem (11.11) in the form of a small perturbation of V, the hrst approximation of (12.201) . 
Since we have, by Lemma 12.41 a small error term R, the equation for the perturbation is a linear one with right hand 
side given by this error term R perturbed by a lower order nonlinear term. The mapping properties of this linear 
operator are fundamental in solving for such a perturbation. Not only this, the nonlinearity must remain small if, say, 
an iterative scheme is produced. An obvious way to write this perturbation is in additive way, say w = W + (j) where the 
remainder term (j> will represent a lower order correction, that is 4> goes to zero as p goes to zero. The nonlinearity term 
N{<f)) produced when substituting this new ansatz w in (11.11) is a polynomial in (p carrying at least quadratic terms. 
We aim at finding solutions for (p small provided that the points C* are suitably chosen. 

In term of a small <j), we can rewrite problem dEH) as a nonlinear perturbation of its linearization jCe{(j>) := (—A)"*(() — 
T{y)(j), namely. 


where 


r £ei(l>) = -R + N{(l>) in He 

(2.30) 

( Bjifi = 0, |j| < m — 1 on ditg, 

T{y) := V{ey)e^''^\ 

(2.31) 

R{y) := i-ArW-T{y), 

(2.32) 

iV(0) :=T(y)(e'^-0-l). 

(2.33) 


Let us observe that near the concentration points Ci the linearized operator Ce is a small nontrivial perturbation of the 
polyharmonic operator while it is essentially this operator in most of the domain. 

Now we intend to solve (12.301) . To do so, we need to analyze the possibility to invert the operator jCe{<t>) in order to 
express the equation as a fixed point problem. It is not expected this operator to be, in general, globally invertible 
because, when regarded in the entire this operator does have kernel: functions ICo and Yij, dehned in (13.21) and 

(lOl) below, with j = l,...,2m, i = 1,..., fc, annihilate it. In suitable spaces, for instance L°°, these functions span 
the entire kernel. In a suitable orthogonal to this kernel, the bilinear form associated to this operator turns out to be 
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uniformly positive definite. Then we are intended to solve a suitably projected version of our problem for which a linear 
theory is in order and, after which, the resolution comes from a direct application of contraction mapping principle. 
The next step will be adjust the points ^ in order to have solutions to the full problem. The latter problem will turn out 
to be equivalent to a variational problem in Our main results will be a consequence of solving this Hnite dimensional 
problem. Moreover, we have the validity of the following estimates 

Lemma 2.6. 

For y £ Qe we have that 

(1) ||r|U = o(i), 

(2) ||A^((/)i) — N((f) 2 )\\t < C'maxi=i ,2 ||0i||oo||</>i — </> 2 ||oo (N is a contraction), 

(3) jjAf((^)||* = 0(11011^) as ||i^|loo 0 (N is almost quadratic). 

Proof 

From Lemma 12.31 we have 

2m 

V{ey)e-^y^ = + 0(.|y - ?'!)) 

that is 

T{y) = a2m{2m - l)!To(y)(l + 0(e\y - ^'|)) 
where To{y) ;= X]i=i ■ This fact gives the validity of the first claim. 

For the general problem (—= V{x)f(u) in a domain it is possible to rewrite (12.3211 and (12.3311 in the following 
way, respectively: 


R{y) = {-Arw + V{x)fiW) 

= Vix)[fiW + <t>)-f{W)-f'{W)),] 


where W is the first approximation of the predicted solution and f' is a functional derivative in the Frechet sense. 
Hence, from this, we gain that jC{(j)) = (—+ V{x)f'{W)(j} is the related linearized operator. 

We calculate 


- N{ct,2)\ = Vix)\fiW + <fi)-f{W + ),2)- f'{W){<l>l-)>2)\ 

= Vix) 


= Vix) 


f f'iW + t),i + il-t))) 2 )dt- [ f'iW)i)>i-(j, 2 )dt 

Jo Jo 

[ (01-</>2)dt+ f /"(IT + st(/)l + s(l - t)02) ■ 

^0 Jo 


it(l>l + (1 + t)cj) 2 ))ds 


and by Lagrange’s theorem we have 

< sup|/"||(/>i - (/ 2 I [ \t(t)i - (1 - t)(j) 2 \dt < sup|/"||(()i - (j) 2 \ max(|(/i|; 1021 ). 

Jo 

This fact gives the validity of the second claim. 

Now, using this second claim with 0i = 0 and 02 = 0 we, easily, get the third. 

□ 


Remark 2.7. 

(1) For the previous analysis, using 2., 3. of Lemma 12.61 we can prove that also N — R is a. contraction. In fact, 

||(iV(0i) -R)- (iV(02) - R)\U = ||iV(0i) - iV(02)||* < Cmax ||0i||oo||0i - 02||co 

1 = 1,2 

and, recalling that R is small by Lemma [2.41 

||(iV(0) - R)\U < ||iV(0)||* + ||R||* < CII 0 IIL + 0(8). 

(2) We note that far from the points (,i (i.e. on most of the domain Q.) T{y) = hence Cei4‘) is a small 

perturbation of (—A)"* away for the concentration points while near the points T{y) — {y)), 

hence Te(0) is a nontrivial perturbation of (—A)"* thus Te(0) is approximately a superposition of the linear 
operators £i(0) = (—A)"‘0 — e™‘0- 
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3. Analysis of the linearized operator 

A main step in solving (12.3011 for small </> under a suitable choice of the point is that of a solvability theory for the 
2m-order linear operator under suitable orthogonality conditions. We shall devote this section to prove the bounded 
invertibility of Le- We consider in 

:= (-Ar^-r(y)0 

for function tj) defined on Qe and where T{y) was introduced in 112.3111 . 

Unlike the 2m- Laplacian (—A)"*, the operator Ce as an approximate kernel which in principle prevents any form of 
bounded invertibility. In fact, centering the system of coordinates at, say, by setting z = y — one can easily see 
that formally Zlg as e —>■ 0 can be approximately regarded as a superposition of the linear operators in 

^ ^ ..‘2m 

r,( 0 ) := (-A)"*^ - ^ e -^0 = (-A)"*^ - a2m{2m - 1)! ^ (3.I) 

i=l i=l + 1^1 > 

namely, Liouville equation (—A)”*in = e™ linearized around the standard bubble 

Wiiz) = a2m{2m - 1)! ^ log (^2 ■ 


Thus the key point to develop a satisfactory solvability theory for the operator jCg is the non-degeneracy of the solutions 
Wi up to the natural invariances of the Liouville equation under translations and, due to the presence of the critical 
exponent, also under dilations, that is if 2; is a solution, necessarily also the function x 1—>■ is a solution, for 

any y > 0 and for some suitable constant a. 

Non-degeneracy is an important ingredient in the construction of solutions to problems involving small parameters 
and concentration phenomena in which, after a suitable blowing-up around a concentration point, one sees a limiting 
equation. This property, in general, is used to build solutions with multiple concentration points. We are interested in 
the classification of bounded solutions of = 0 in R^™'. Some bounded solutions are easy to find. For example, we 

can define the functions 

where r = \z\ and, clearly, Ci(Yio) = 0 and this reflect the fact that the equation is invariant under the group of 
dilations r —>■ —u(r-) -|- 2m log r and for all j = 1,..., 2m and i ~ 1 ,..., k 

Yij{z) =-d^^Vi = Am 2^1,12 ’ (3-3) 

k-i I \z\ 

clearly CiiYij) = 0 and this reflect the fact that the equation is invariant under the group of translations a —>■ «(• -|- a). 
In the next Lemma, following very close the geometric point of view in Lemma 3.1), we prove that any solutions of 
the linearized operator Ci calculated along the standard bubbles are nondegenerate, that is it has not trivial kernel in 
certain space, for instance L°° , and thus, by invariances of the Liouville equation under a large groups of symmetries, jCi 
has bounded kernel. This is equivalent to the fact that the actually L°°-kernel of the linearized operator ker jCir\L°° (R^™') 
is spanned by the bounded solutions (13.21) and (13.31) which naturally belong to this space, that is the only bounded 
solution of = 0 in all R^™' are linear combinations of those functions. In the proof, we use some well-known 

fact in spectral theory and the classical idea to prove nondegeneracy comparing precisely the dimension of the bubble 
solutions manifold with respect to the dimension of the kernel of the linearized operator calculated along the standard 
bubble. As a byproduct, we have that ker£i is equal to the eigenspace of Pgg2m associated to the eigenvalue tm and 
that dim(ker/li) = 2m -|- 1. 


Lemma 3.1. (i°° non-degeneracy) 

Any bounded solutions of = 0 in all R^"* are linear combinations of (13.21) and (13.31) . 

Proof 

We consider on R^™' the Euclidean metric gs = \dz\'^ and the spherical metric on 

induced by the inverse of the stereographic projection from the sphere onto the whole R^"* with respect to the 
North pole, namely 

: 2 = (zj) G R^"* \ {North pole} C R^™'+^ 

given by the formulas 

l<i<2m 


^(2) := 


T+jfl^ j = ‘2m + l. 


We remark that tt ^ is a conformal diffeomorphism, more precisely the pullback of to R^"* satisfies 


(t ) Ps2m — 


(1 + kr 
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According to [6l|8], we have Pg^sm = + k{2m — k — 1)) where Pg^sm i® ^he G'JMS'-operator on and 

— m is the Laplace-Beltrami operator on and when the manifold is the Euclidean space, the G'JMS'-operator is 
simply given by PR2m = (—A)"*. Since under the following conformal change of metric p = the GJMS-operator 

transforms according to Pg = e“^™''^Pg(/) we obtain that I j P9s2m = Pse- 


In particular, if </> = ^ovr : (R^"*, (;_b) —>■ R (with slight abuse of notation we identify (f> with (/iott) is a bounded solution 
of Ci{4>) = 0 then 0 = 0 o tt : (R^"*, g^ 2 m) —>■ R is a bounded solution of 


P 9 s 2 m ^ in \ {North pole} 


(3.4) 


away from the North pole and where tm ’■= = 2m(m — l)!(2m — 1)!. For instance, t 2 = 2, t 4 = 24. 

Since cf) is assumed to be bounded then the isolated singularity at the North pole is removable and hence (E3 holds 
on all S^"". 

We now perform the eigenfunction decomposition of 4> in terms of the eigendata of the Laplacian on S^"*. We decompose 


0 = 

k >0 


where (f>k belongs to the fc-th eigenspace of Ag 2 m, namely (fik satisfies the following eigenvalue’s problem —^g^ 2 m 4 ‘k = 
Xk 4 ‘k- It is well known that the spectrum of — Ag^j^ is discrete and represented by 

<^(-^ 9 ^ 2 ™) = A > 0} = {k{k + 2 m- l)/k > 0} 

with multiplicity ruk = ^ gee f|25|. Appendix C). By (13.41) we have 


m—1 


n 


(—Afe + k{ 2 m — k — 1)) — tm 


(pk = 0 , 


that is YYk^o^~‘^k'^ — tm)<t>k = 0. Hence (j>k =0 for all fc > 0 except, eventually, those for which YYk=o{~‘^k'^ — tm) = 0. 
This implies that (j) : ^ R is a combination of the eigenfunctions associated to fc = 1 that are given by ipj (y) = yj , 

j = 1,..., 2m + 1, where y — (yj) G Being the sphere parameterized by 7r“^ we may write y = 7r“^(2:). Then, the 
functions 2 my 3 j precisely correspond to the functions Yij while the function 2 mifi 2 m+i corresponds to the function Yio- 

□ 


4. Projected linear theory for Ce onto kernel 
We define for j = 0,..., 2m and * = 1,..., k, 

^ijiy) •= 4^0 (^)' 

Additionally, let us consider Rq a large but fixed number and x : R ^ R) Xiiv) = x(Ti)i where Vi := \z\, a smooth and 
radial cut-off function with x := 1 if Ci < J?o, namely in Bho (0) and x := 0 if ri > i?o +1, namely in the complementary 
B«o+i(0)^ 

Given h G L°°{Qe), we consider the linear problem of finding a function <() : > R, <() = (fi{y), such that for certain 

scalars dj one has 

( = h + CijXiZij in Xle, 

< Bj(j} = 0,\j\ < m — I ondQe, (4.1) 

I Ifie = 0 for all j = 1,..., 2m, i = 1,..., k. 

The orthogonality condition in (14.11) are only taken with respect to the elements of the approximate kernel due to 
translations. 

The main goal of this section is the bounded solvability of dm. Before this, we will establish a priori estimates for 
this problem. To this end we shall conveniently introduce an adapted norm in fig. Given (/>: fig —> R and a G N*’ a 
multi-index of order k and length |q| we define 

k k 

IAII« :=EIAIIc-,<^(..<2)+E E (4.2) 

i=l i=l |Q:|<2m —1 

We remark that the interior portion of the norm defined above controls the function <() in a neighborhood of the origin 
while the exterior portion the decay at infinity of (j). With this definition at hand, we prove the main result of this 
section: 

Proposition 4.1. 

There exist positive constants eo and C such that for any h G with ||h||* < oo, and any ^ G A4, there is a 

unique solution (j) = Q{h) to problem (14.11) for all e sufficiently small, say e G (0, eo), which defines a linear operator of 
h. Besides, we have the a priori estimate 


lAII** < C|loge|||h||*. 


(4.3) 
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Before proceeding with the proof we remember a useful lemma on the removable singularities for poly harmonic function. 
Denote by B and Bq the n-unit ball and the n-punctured unit ball B \ {0}, respectively. We present the result in ( |29|. 
Theorem 3.2) specified to the n = 2m case. 

Theorem 4.2. 

Suppose n = 2m and u a m-harmonic function defined on Bq. Then, the following are equivalent: 

(1) u can be extended to a m-harmonic function on B; 

(2) lim u{x) exists and is finite; 

cc —>-0 

(3) u is bounded near the origin. 


The proof of Proposition 14.11 will be split into a series of lemmas which we state and prove next. The first step is to 
obtain a priori estimates for the problem 

( = h in De, 

< Bj(j) = 0,\j\ < m — 1 ondQe, (4.4) 

I fcie XiZij4> = 0 for all j = 0,..., 2m, i = 1,..., k 

which involves more orthogonality conditions than those in (HU- Notice that in the case m = 1, independently from 
the nonlinearity, a key step in order to prove such result is the fact that the operator Ce satisfies maximum principle 
in Tie outside large balls, see ([23], Lemma 3.1) and |27| for the exponential-type nonlinearity or the surveys [49] and 
references therein for other nonlinearity issue such as the Brezis-Nirenberg Problem or the Coron’s Problem. For our 
Ce this is not more true and we need a different approach. We have the following estimate 

Lemma 4.3. 

There exist positive constants eo and C such that for any solution (j> of problem (14.411 with h G L°°{Tle), II ^||* < oo 
and with f, £ M, then for all e € (0,£o) 

11011** <C||/i||*. (4.5) 


Proof 

We carry out the proof by means of a contradiction argument. If the above fact were false, then, suppose there exist a 
sequence —>■ 0, a sequence of points £ M with i = 1,..., k, a sequence of functions h„ with ||h„||* —>■ 0 as 

n —>■ 00 and a sequence of associated solutions fin with ||0„||** = 1 of the following problem 

{ Cen {.fin) ~ hn iu Tle^. , 

Bjfin = 0,1 j| < m - 1 on dQen, (4.6) 

fcie XiZijfin=0 for all j = 0,..., 2m, i = 1,..., fc, 

where Ce^ifin) := (-A)"*(^„ - V{eny)e'^^^^fin- 

We observe that ||/in||* —>■ 0 as n —>• oo because if (14.51) doesn’t hold then ||((>n||** > C||/irt||*, normalizing this last 
inequality dividing by ||(^n||** we have done. We will show that, under those hypothesis, we necessarily have fin ^ 0 in 
Qe and, then, a contradiction arises. 

Let us set fin := fin{Sn^x), X £ fl. First, we prove an intermediate claim: with the above hypothesis we have that 
fin goes to zero as n —>■ oo in C'™'^^’“-sense uniformly over compact subsets of D \ {^*}. In particular, for any J' > 0 
sufficiently small we have 

\\fin llLoo(n\ufc^,^ (?„) ^ 0 as n —>■ oo. 

Now, fin solves the equation in : 

{-A)'^fin{x) = en^”"{-A)""finien^x) = En^"" [-Tn{efi^ x)fin{en^ x) -|- hn{£n^ x)] 

= £-^^[V{x)e'^^^-"^'^fin{x) - hn{£fi^x)] = [V {x)e^ ^^^n{x) - hn{£n^ x)]. 

For any <5' > 0 sufficiently small, we recall that in \y — (^")^| > £n^5' where {£fi)' '■= we have £^Tn{y)fin{y) = 

0{£^) because in this region Tn{y) = 0{£^), that is small, and fin is uniformly bounded by contradiction hypothesis. 
Setting X = £„?/ and recalling that W{y) = U{x) -I- 2mlogpe„, then also e^V{x)e^^^'*fin{x) = 0{e^). Thus, 
{-A)^fin{x) = 0(4"^) + e-2"‘h„(£-4). 

At this point, we claim that 0{e^) -£ £n^^hn{£n^x) = o(l) as n —>■ oo uniformly in D := fl \ IJiLi Bs'{^n). In fact, 
if \x — ^„| > S' then the dominant part in the definition of || • ||, is, as noted before, the second, thus, £))'^’”|h„(j/)| < 
c||hn(2:)||, and by contradiction hypothesis we easily get the desired result. We have directly checked that for any <5' > 0 
small fin solves the problem 

f {—A)"'fin = o(l) uniformly in fl, 

\ Bjfin = 0, |ji| < m — 1 on 9fl, 



SINGULAR LIMITS IN HIGHER ORDER LIOVILLE-TYPE EQUATIONS 


17 


together with the sharp estimates: 

< 1 and < Cis j = 1, 1. (4.7) 

Differently from m, Lemma 3.1) we are also able to compute the constant Cs'- 

Fix an index 0 < q < 2m — 1, since by contradiction hypothesis ||(()n||** = 1 then the interior part of the norm || • ||** 
vanish while for the exterior part we have that 

||zrDV„(2)| < 1, V2 := G ne„ \ 52. (4.8) 

Since 4in{z) = 4>n{x) then D'^(f>n{z) = <j)n{x). Hence by (14.811 we get 

||a:|'^5‘^<^„(a:)| < 1. 

If g = 0 this means that ||0n(a:)||i=o(n\s2e„) — 1 since Eniyie^ \ B 2 ) = 0. \ 52e„. Now, if g > 0 and x a fl \ 5^/(0) then 
we have (S')'^\D'^(f>n{x)\ < \\x\'^D'^(f>n{x)\ < 1. If we replace q with 2q we have Finally, 

since this is an estimate for each q we have to sum over q and thus we introduce a constant Cq, dependent only from 
q, in order to have the desired estimate with c^/ := jjt^- Notice that, since M is compact, we may always pass to a 
subsequence, that with slight abuse of notation we continue to call := {^h), such that ^ C* G as n —^ 00. Then 
using by estimates (T71l and by m, Theorem 2.6) there exists a subsequence still denoted by cj)n such that (jjn 4>o as 
n —>■ 00 in C^™'“^’“-sense over compact subsets of fl \ {^*} because > fl \ {^*} as n —^ 00. Furthermore, ^0 solves 

r (-A)"‘,^o = 0 inL!\{r}, 

\ Bj(l>o = 0, \ j\ < m — 1 on 


By the first estimate in follows that 
thus, solves 


< 1- By Theorem 14.21 then the singularity is removable and, 

(-A)™-00 = 0 in n, 

5j0o = 0 , IjI < m — 1 on dil. 

If m = 1, m > 1 with Navier b.c. or m > 1 even and Dirichlet b.c. integrating by parts we have directly that 0o = 0. 
If 1 < m is odd and with Dirichlet b.c. in order to conclude we need also to use, after an integration by parts, the 
well-known decomposition of the Laplacian on dfl valid for any smooth function u, namely 

Au = dlu + HendvU + AtU on 9D, 

where Hgn 7^ 0 denotes the mean curvature at the boundary and At denotes the tangential Laplacian, see ([55], 
pag.62). For instance, if m = 3 using by the Dirichlet b.c. the above formula reduces to Au = d^u on dQ,. Finally, up 
to a subsequence, 0„ —>■ 0 as n —>■ 00 in C^"‘“^’“-sense over compact subsets of D \ {C*}. In particular, 

^ (e"^)'“'|D“0„(y) 0, uniformly in \y - (5„)'| > 

|a|<2m —1 ^ 

for any > 0 and 2 = 1, k. We obtain thus that 

k 


E E IIE'5“ 

i = l \a\<2m—l 




0, 


(4.9) 


for any S' > 0. In conclusion, the exterior part of the norm ||0„||** goes to zero, thus the claim follows. 

Let us consider now a smooth radial cut-off function 77 with rj{s) = 1 if s < ^,ri{s) = 0 if s > 1, and dehne 


such that 

Notice that 
where 


0n(l/) := v(^^\y- (Cnl^'('n(y), 

supp0„ c B^{^'„). 

-Cen (0n) = vhn + F{r], 0„), 


F{f,g) (-A)’"(/5) + C^ 




d'^f d'^g 
dyidyj dyidyj ' 


Thus we get (—A)™’0„ = Tn{y)^n + rjhn + F(y, (jin) in B^ (^(,) and, since 0n G Co°{B^ (^n)), in particular we have 

En fin 

that Bj(f>n = 0 on dQ. Since ||0n||** = 1, using (14.911 and the outer estimate in ([19], Lemma 3.2) we have that there 
exists, up to a subsequence if necessary, an index i = 1,..., fc such that 


lim inf | 

n—^oo 


l|n“(ri<i?o) > « > 0. 


(4.10) 
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Let us set ipn , i { z ) := 4 > n{z + where the index i is such that sup|^_j/ |0„| > a > 0. Without loss of generality, 
we may always assume that the index i is the same for all n . Notice that satisfies 

{-A)”^1pn,i - T„{z + i'n)tpn,i = hn{z + iu Q.n := - {Cn}- 

Since Bj(j)n are bounded uniformly then as n —>■ oo over compact subsets of we have 

Tn{z+a = ^^^l^^ii+o{i)), \h„{z+a\<c\\h„\u. 

Hence, standard elliptic estimates allow us to assume that, as n —^ oo, ipn,i converges uniformly over compact subsets 
of to a bounded, non-zero solution ip of 


(-Arv' = 


a2m(2m - l)!/rf 


-Ip in ' 


This implies that, by Lemma [3.11 -i/i is a linear combination of the functions Yij defined in (13.211 and (13.311 . Orthogonality 
conditions over ipn,i pass to the limit thanks to ||oo < 1 and dominated convergence theorem. Thus, this implies 
that tp(y) = 0, a contradiction with (14.1011 . 


Now we will deal with (14.41) lifting the orthogonality constraints Xi^io<P = 0, i = 1, namely 

( Ce{(p) = h in He, 

< Bjtp = 0, |j| < m — 1 on 9He, 

I /oe = 0 for all j = 1,..., 2m, i = 1,..., k. 

We have the following a priori estimates for this problem. 

Lemma 4.4. 

There exist positive constants eo and C such that for any solution (p of problem (14.1111 with h G L°°{Qe), 
and with f £ A4, then for all e £ (0, £o) 

||<^||** < C\ logel 


□ 


(4.11) 


* < 00 
(4.12) 


Proof 

Let i? > 7?o + 1 be a large, fixed number and Rq is the radius of the ball where the cut-off function x, defined in the 
previous Lemma, is supported. Let us consider 


Zio{y) ■- Zoi{y) - 1 -|- aioG{ey,^i), 


(4.13) 


where aw := — 4mlog(eR)) Notice that, by (11.911 . aw = ^ small enough then, by lESl) and 

Lagrange Theorem we have 


(4.14) 


^io(y) = ^io(y) + awl G{ey,ii) - + 4mlog{eR) 


— Zw{y) + 


logs I 


0{eri) -I- 4m log 


R 


and Zw{y) = 0(1). Now, we consider radial smooth cut-off functions 771 and 772 with the following properties: 

0 < 771 < 1, 771 = 1 in Br{0), ryi = 0 in (i3H+i(0))‘', 

and 

0 < 772 < 1, 772 = 1 in Bso (0), 772 = 0 in (Bsg (0))'". 

3e 2e 

Without lack of generality, we may assume Ssq ( 0) C Q. Set 

Vii{y) ■- Vi{ri), 7712 ( 7 /) := 772 ( 7 - 1 ), 

and define the test function 

Zw '■= ViiZw + (1 — r]ii)rji2Zw- 

Notice that suppZio C H and are all disjoint. Furthermore, intuitively, Zw resembles the eigenfunction of the operator 
Ce in R^™- with respect to the dilation property of C-e- Observe the Zw's, behavior through Tig-. 

Zw in Ho {ri < r}, 

rjii {Zw — Zw) + Zw in Hi {i? < < i? -|- 1}, 


Zw = > 


Zw 
7712 Zw 
0 


in H2 ;— {i? -|- 1 < ri < ^}, 
in Ha := {fj < < |j}. 


otherwise. 
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Let 0 be a solution to din . The main idea of the proof is to modify 0 so that the extra orthogonality conditions with 
respect to Zio’s hold and try to use Lemma 14.31 Set 


k 


0 :— 0 -h ^ ^ diZiQ. 
i=l 

(4.15) 

Our goal is to adjust the constants di so that 


[ XiZij^ = 0, Vj = 0,2m; i = 1, k. 

J Qg 

(4.16) 

Then, 


k 

jCe{4>) = h + ^ diC,{Zio) in Tl,. 

(4.17) 


i=l 


If (14.1611 holds, Lemma allows us to conclude 

||<^||** <c|||h||*+^|d,|||£e(^io)||*|. (4.18) 

i=l ^ 

Estimate din is a direct consequence of the following claim. 


Claim: The constant di are well defined, 

\di\ <C\ log£|||/i||* and ||£e(^oi)ll* < ,. ^ , , Vi = 1,..., k. (4.19) 

|log£| 

In fact, using by this claim and the fact that ||Zio||** < C we obtain din as desired. Now, let us prove the claim. 
First we find di. From (14.1511 . orthogonality conditions (14.161) and the fact that suppxj77ifc = 0 and suppxj772fc = 0 if 
i 7^ j are also satisfied for thanks to the fact that R > Ro + 1 we can choose 


di = 


In, XiZio<l> 
In, Xi\Zio\^’ 


Vi = 1,..., k 


and, then, di is well defined. In order to prove the second inequality in (14.1911 we have to compute jCs{Zoi) in fl;, I = 
0, ...,3. Proceeding exactly as in m, Lemma 3.3, Claim 1) we may obtain the desired estimate. Finally, known this, 
we may prove the first inequality in (14.191) . 

Testing equation (14.171) against Zio and the above estimate, we get 



£.e(Zio)Zio 


hZio + 



<C\\h\\,+CU\\^\\CeiZio)\U. 

Since ||Zio)||* = 0(1) Vi and ||(?i||oo < ||'('||»» using by relations (14.181) and the second inequality in (14.191) we have 


\di\ 


I 

Jq. 


C£{ZiQ)Zio 


<c\\h\u+cJ2- 


\di\ 


logel' 


(4.20) 


It only remains to estimate the integral term in the left side. Proceeding exactly as in m, Lemma 3.3, Claim 2) we 
may obtain the following claim. 


Claim: If R is large enough, then 


£e{Zio)ZiCi 


C 


|loge| 


( 1 + 0 ( 1 )). 


(4.21) 


where C is a positive constant independent of e and R. 


At this point, we may replace (14.211) in (I4.2UD in order to get the desired bounds of di. 


□ 


Proof of Provosition \4:.l\ 

First we establish the validity of the a priori estimate (|43J for solution if of (14.11) with h G L°°{Q,e) and ||h||* < oo. 
Lemma 14.41 implies 

k 2m 

||<^IU*<qiog£|{||h|U + ^^|c,,|||x,z,,|l*}. 

i=l j = l 


(4.22) 
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Since £ C then it sufficient to estimate the values of the constant dj. In order to do this, we multiply the 

first equation in 631) by ZijTju, with rii 2 as in Lemma 14.41 and integrating by parts to find 


/ 

Ja, 


I 


— / hZijT^i2 ~\~ Cij I jZijj Tji2- 


I 

Jfl. 


It is easy to see that hZijrji^ = 0(||/i||*), = C > 0 and 


(4.23) 


/ 

Ja, 


C,eZijTji 24 > 


< Ce\ log£|||0||co < Ce\ logE 


see also [I1123111I]. Using the above estimates in (14.2311 we have 


\cij\ < C{e\ log£111011** + ||/i||*}, 


and then we get 


'{(1 + 


£| l0g£|' 


+ £|l0g£|^y]| 


l,n 


Thus, |cij| < C||/i||* and putting this estimate in (14.2211 we conclude. 

Now, we prove the solvability assertion. Notice that Problem (14.1111 expressed in a weak form is equivalent to that of 
finding a 0 G 'H, such that 

(0,0)_ff= [ {h + T(j))tp, yip eK, 


where we consider the Hilbert space 


H ~ {(p '■ Bj(p = 0 on dQ,e and 



XiZij(p = 0, Vi = 


1, ...,fc;j = l,...,2m}. 


endowed with the usual inner product {(p,ip)H. With the aid of Riesz’s representation Theorem, this equation can be 
rewritten in 'H in the operator form (p = K{T(p + h), where R' is a compact operator in H. Then Fredholm’s alternative 
guarantees unique solvability for any h provided that the homogeneous equation <p = K(Tip) has only zero solution in 
H. This last equation is equivalent to (l4dTll with h = 0. Thus existence of a unique solution follows from the a priori 
estimate (14.1211 . This finishes the proof. 


□ 


Remark 4.5. 

(1) This result implies that the unique solution 0 = Q{h) of (14.11) defines a continuous linear map from the Banach 

space C* of all functions h € with ||/i||* < +oo, into with norm uniformly bounded in 

£. 

(2) The operator Q is differentiable with respect to the variables . In fact, computations similar to those in ([23], 
pag.l7) yield the estimate 

||%Q(fe)||** <qioge|"||h||*. (4.24) 

This estimate is of crucial importance in the arguments to come. 


5. The intermediate nonlinear problem 


Rather than solve Problem (I2.3UI) directly we shall consider the intermediate nonlinear problem 
( ^ei<P) =-R +N{(p)+ Y^‘^.'PP^Y^'p^j^CijXiZij in He, 

< i3j0 = 0, |j| < m — 1 on 5He, (5.1) 

I In, XiZijfp = 0 for all j = 1,..., 2m, i = 1,..., k. 


Assuming that conditions in Proposition 14.11 hold, we are able to prove the following 

Lemma 5.1. 


Let ^ G A4. Then, there exists positive constants eo and C such that for all e < eo the nonlinear Problem 631 has 
a unique solution 0 which satisfies 

II0II** < Ce|log£|. 

Moreover, if we consider the map ff £ M ^ (p £ C^"‘’“(He), the derivative D^Kp exists and defines a continuous map 
of f'. Besides, there exists a positive constant C such that 


D^,(p\\„ < Ce\ log£|^. 


Proof 

In terms of the operator Q defined in Proposition 14.II (15.11) has the following fixed point representation 


f> = B{cP) = Q{N{cP)-R). 
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Let us consider the region 

||</.||„ <£|loge|}. 

From Proposition O we have 

||S(0)||**<qioge|{|iV(0)|U + ||7?|U}, 
for arbitrary </>. Hence, by Lemma 12.41 and Lemma ESI we have V(/>, (/>i, 02 £ 

(1) 11^(0)11** < C'el log e|; 

(2) ||H(0i)-H( 02)||** <Ce|logep||0i-02||**. 

Then it follows that for all e sufficiently small H is a contraction mapping of T, and therefore, by using the Implicit 
Function Theorem, a unique fixed point of B exists in this region. The IFT guarantees regularity of the map in 
Follows exactly the same proof included in m, Lemma 4.2) we may show the derivative estimate. 

□ 


6. Variational reduction 

After Problem eH) has been solved, we will find solution to the full Problem (12301) . or equivalently (ini), if we 
manage to adjust the fc-uple in such a way that 

co (C0 = 0i for alH, j. (6.1) 

A nice feature of this system of equations is that it turns out to be equivalent to finding critical points of a functional of 
^ which is close, in appropriate sense, to the energy of the first approximation W. Notice that problem (16.111 is indeed 
variational. In fact, to see that let us consider the energy functional Jp associated to Problem (11.11) . namely (12.251) . We 
define the function for ^ a M 

5e(O = ^p[l^(O + 0e], (6.2) 

where U = U{^) is our approximate solution from (11.71) and 0^ := 0(-, -), x € and 0 the unique solution to Problem 
(15.11) predicted by Lemma l5. II Under the assumptions of Lemma EH we obtain, in the following result, that critical 
points of correspond to solutions of (EH for small e and, furthermore, the closeness of to Jp[U(^)], for which we 
know the asymptotic estimate (12.261) . 

Lemma 6.1. 

The functional : A4 —>■ K is o/ class . Moreover, for all positive constant e small enough, if D^iFeif.) ~ 0 then 
f satisfies EH- Besides, for ^ G A4 the following expansion holds 

M^,) = JAUm+0e{i), (6.3) 

where \D^^0e\ = o(l), uniformly on £ M os e —^ 0, |jj < m — 1. 

Proof 

We define the functional 

^elw].= ^[ \{-A)^w\'^dy- f V{£y)e^dy. 

^ Jcie JOe 

Let us differentiate 'Se with respect to 5. Notice that, since ./p[U(5) + 0{] = Ie[W{^') + we can differentiate directly 
under the integral sign, so that 

%,)n5e(0 = ii3/e[W + 0](%,)„W + %,)„0) 

= 7Y.Y.I +%,)»<(>)• 

From the results of Section El this expression defines a continuous function of (f, and hence of f. Let us assume that 
= 0. Then, 

/ cuV^u(%,)nW^ + %i)n'(') = 0. / = l,...,2m,n = l,...,fc. 

3 i 

Since, by the derivative estimate in Lemma EH] we have directly 

where o(l) is uniformly small as e —>■ 0. Thus, we get that = 0 implies the validity of the equations in term of 

**-norm 

y + o(l)) = 0. 

j i 

This system is dominant diagonal, thus we get dj = 0 for all i,j. 

The closeness property follows directly from an application of Taylor expansion for Je in Qg and from the estimates in 
Lemma EH See also ([23], Lemma 5.2). This concludes the proof. 

□ 
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7. Proof of Theorem o 


Taking into account Lemma l6.ll a solution to uni) exists if we prove the existence of a critical point of which 
automatically implies that dj = 0 for all i,j. Next, the qualitative properties of the solution found follow from the 
chosen ansatz. Finding critical points of is equivalent to finding critical points of 

de ■= de -‘imbm.k] logs]. (7.1) 

On the other hand, if ^ a M, from Lemma [23] and Lemma IQ we get the existence of universal constants a > 0 and 
/3 such that 

ade+l3 = M0+O{e). (7.2) 

With the same argument in m, Theorem 2) and ([23], Theorem 2), we may prove that, under the assumptions of 
Theorem OSe has a critical point in M for e small enough. By 113 the proof is concluded. 

□ 


Appendix 

In the sequel we calculate the explicit values co and ci used in Lemma [2.51 This values, appeared until now only in the 
case m = 1, 2, are broadly used when we work with elliptic equations with exponential nonlinearity. 


Lemma 7.1. 


/'I'l f d-y _ 7r™-(m-l)! 


(r)\ r log(l+|v|^)dy _ 7i-™(m-l)! 
Jll2m (I+IhI^)^’" 

Proof 


—1)! 
m(2m—1)! 


We start to calculate the following integral /jj 2 m for ot > m, interesting by itself. The formula of integration 

in spherical coordinates of an integrable rotationally symmetric function yields 


L 


dy 


We observe that 


L 


R 2 m (l + |l/P)“ 

da 1 


/ 8 fl{ 0 ,r) (l + r2)“ (l + r2)' 

By the change of variables r = faxiO, 


ru 

Jo \JdB(0,j 


-2muJ2mr 


da 


) (1 + r2)= 


dr. 


r(m)(l + r 2 )° 


rcG 2m —1 

I = 1 (cose)^(“-)-^(sine)= 


^de. 


Recalling the definition of the Beta/unction/J(s, t) := 2 ^(sinS)^* for all (s, t) £ , the set of positive 

rational numbers, we have 


L 


oo ^2m-l ^ 

j^^-^dr = -l3ia-m,m). 


Recalling the relation between Beta and Gamma function, we have JJ” Tr-rr’rradr = 4 and 


thus, 


/ 


dy 


r(s+t) 

7r’"r(a — m) 


Jo (l+r2)o 


r(c) 


/R2m (1 + |i/|2)“ r(a) 

In particular, with a = 2m and by Legendre duplication formula for T: 


we have 


/ 


r(m)r^m+i^ =2^ ^/TvV^m) 
7r™'r(m) 


dy 


2-k^ 


/R2m (l + |yP) 2 - r(2m) 22-r(m + i)’ 

Using the classical equality r(m+ i) = , the fact that ( 2 m— 1)!! = and by the definition of the factorial 

of a nonnegative integer the proof of the first point is concluded. 

Moreover, integrating by parts, we note that 


log(l + \y\'^)dy _ 2 muj 2 r, 


Vm (1 + |i/|2)2" 

and the proof is concluded. 


L 


dy 


m V-n (1 + | 2 /P) 2 ’" 


= 2lJ2t 


/. 


dy 


R2m (1 + 1^12)2"* 


□ 
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